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m ■ 

We consider the evolution of small rotational perturbations, with azimuthal symmetry, of the 
brane-world cosmological models. The equations describing the temporal, radial, and angular de- 
£Sj ■ pendence of the perturbations are derived by taking into account the effects of both scalar and tensor 

parts of the dark energy term on the brane. The time decay of the initial rotation is investigated for 
O ' several types of equation of state of the ultra-high density cosmological matter. For an expanding 

Universe, rotation always decays in the case of the perfect dragging, for which the angular velocity 
of the matter on the brane equals the rotation metric tensor. For non-perfect dragging, the behavior 
of the rotation is strongly equation of state dependent. For some classes of dense matter, like the 
stiff causal or the Chaplygin gas, the angular velocity of the matter on the brane is an increasing 
function of time. For other types of the ultra-dense matter, like the Hagedorn fluid, rotation is 
smoothed out by the expansion of the Universe. Therefore the study of dynamics of rotational 
perturbations of brane world models, as well as in general relativity, could provide some insights on 
. the physical properties and equation of state of the cosmological fluid filling the very early Universe. 

o 

oo 
O 
<N 
O 

I. INTRODUCTION 

Most of the astronomical objects in the Universe (planets, stars or galaxies) have some form of rotation (differential 
or uniform). Hence the possibility that the Universe itself could be rotating has attracted a lot of attention. But 
even that observational evidences of cosmological rotation have been reported Ql d, HI, Q , they are still subject of 
controversy. 

From the analysis of microwave background anisotropy Collins and Hawking 5], and Barrow, Juszkiewicz, and 
Sonoda have found some very tight limits of the cosmological vorticity, T t, s > 3 x 10 5 Th, where T t, s is the actual 
rotation period of our Universe and Th = (1 ~ 2) X 10 10 years is the Hubble time, corresponding to an angular 
velocity of the order of 10 -13 rad/sec. Therefore our present day Universe is rotating very slowly, if at all. However, 
the existence of such a small rotation, when extrapolated to the early stages of the Universe, could have played a 
major role in the dynamics of the early Universe, and possibly as well in the processes involving galaxy formation. 

From a theoretical point of view, Godel gave in 1949 his famous example of a rotating cosmological solution to 
the Einstein gravitational field equations. Godel also discussed the possibility of a cosmic explanation of the galactic 
rotation 0. This rotating solution has attracted considerable interest because the corresponding Universes possess 
the property of closed time-like curves. 

The investigation of rotating and rotating-expanding Universes generated a large amount of literature in the field 
of general relativity, the combination of rotation with expansion in realistic cosmological models being one of the 
most difficult tasks in cosmology (see 8] for a review of the expansion-rotation problem in general relativity). Hence 
rotating solutions of the gravitational field equations cannot be excluded a priori. But this raises the question of why 
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the Universe rotates so slowly. This problem can also be naturally solved in the framework of the inflationary model. 
Ellis and Olive [j|, and Gr0n and Soleng 01 pointed out that if the Universe came into being as a mini-universe of 
Planck dimensions and went directly into an inflationary epoch driven by a scalar field with a flat potential, due to 
the non-rotation of the false vacuum and the exponential expansion during inflation the cosmic vorticity has decayed 
by a factor of about 10~ 145 . 

The rotational perturbations of a spatially homogeneous and isotropic Universe in terms of a small variation of 
curvature have been investigated by Hawking . He found that for pressureless dust and radiation the perturbations 
die away. Later Miyazaki [12| studied the perturbations by a spherical shell in a closed homogeneous Universe in the 
framework of the Brans-Dicke theory. The possibility of incorporating a slowly rotating Universe into the framework of 
Friedmann- Robertson- Walker (FEW) type metrics has been considered by Bayin and Cooperstock 0] , who obtained 
the restrictions imposed by the field equations on the matter angular velocity. They also shown that uniform rotation 
is incompatible with the dust filled (zero pressure) and with the radiation dominated Universe. Moreover, Bayin [l4j 
has obtained the solutions of the field equations for a special class of non-separable rotation functions of the matter 
distribution. The investigation of the first order rotational perturbations of flat FRW type Universes proved to be 
useful in the study of string cosmological models with dilaton and axion fields 0] . The form of the rotation equation 
imposes strong constraints on the form of the dilaton field potential U, restricting the allowed forms to two: the trivial 
case U = and the exponential type potential. 

Recently, Randall and Sundrum [To . Il7| have pointed out that a scenario with an infinite fifth dimension in the 
presence of a brane can generate a theory of gravity, which mimics purely four-dimensional gravity, both with respect 
to the classical gravitational potential and with respect to gravitational radiation. The gravitational self-couplings are 
not significantly modified in this model. This result has been obtained from the study of a single 3-brane embedded 
in five dimensions, with the 5D metric given by ds 2 = e~f( y ^r] l j lV dx fl dx L ' + dy 2 , which can produce a large hierarchy 
between the scale of particle physics and gravity, due to the appearance of the warp factor. Even if the fifth dimension 
is uncompactified, standard 4D gravity is reproduced on the brane. This model allows the presence of large or even 
infinite non-compact extra dimensions. Our brane is identified to a domain wall in a 5-dimensional anti-de Sitter 
space-time. 

The Randall-Sundrum (RS) model was inspired by superstring theory. The ten-dimensional Eg x Eg heterotic string 
theory, which contains the standard model of elementary particle, could be a promising candidate for the description of 
the real Universe. This theory is connected with an eleven-dimensional theory compactified on the orbifold i? 10 x S 1 /Z2 
[l8j| . The static RS solution has been extended to time-dependent solutions and their cosmological properties have 
been extensively studied [H |H M, El H |M |H (for a review of dynamics and geometry of brane universes see 
0). 

The effective gravitational field equations on the brane world, in which all the matter forces except gravity are 
confined on the 3-brane, in a 5-dimensional space-time with ^-symmetry have been obtained, by using an elegant 
geometric approach, by Shiromizu, Maeda and Sasaki |2^,|2^. The correct signature for gravity is provided by the 
brane with positive tension. If the bulk space-time is exactly anti-de Sitter, generically the matter on the brane is 
required to be spatially homogeneous. The electric part of the 5-dimensional Weyl tensor Ejj gives the leading order 
corrections to the conventional Einstein equations on the brane. The effect of the dilaton field in the bulk can also 
be taken into account in this approach |30l |. 

In brane- world model, the behavior of an anisotropic Bianchi type I cosmology in the presence of inflationary scalar 
fields has been considered by Maartens, Sahni and Saini [3j. By using dynamical systems techniques, the behavior 
of the FRW, Bianchi type I and V cosmological models in the RS brane world scenario, with matter on the brane 
obeying a barotropic equation of state, has been studied by Campos and Sopuerta j3l|3i|. The general exact solution 
of the field equations for an anisotropic brane with Bianchi type I and V geometry, with perfect fluid and scalar 
fields as matter sources, has been found in [34|]. In spatially homogeneous brane world cosmological models the initial 
singularity is isotropic, and hence the initial conditions problem is solved [351. Consequently, these models do not 
exhibit Mixmaster or chaotic-like behavior close to the initial singularity [36| . 

Realistic brane-world cosmological models require the consideration of more general matter sources to describe 
the evolution and dynamics of the very early Universe. The effects of the bulk viscosity of the matter on the brane 
have been analyzed in |37| . Limits on the initial anisotropy induced by the 5-dimensional Kaluza-Klein graviton 
stresses by using the CMB anisotropies have been obtained by Barrow and Maartens l3q|. Anisotropic Bianchi type 
I brane- worlds with a pure magnetic field and a perfect fluid have also been analyzed [39| . The effect of the bulk 
viscosity of the cosmological matter on the cosmological evolution on the brane for a Bianchi type I brane geometry 
was considered in |4Cj . 

The simplest way to investigate if brane world cosmologies are consistent with the observations is to investigate 
the behavior of the perturbations in the model. Perturbations on the brane are associated with perturbations in 
the geometry of the bulk space-time. The linearized perturbation equations in the generalized RS model have been 
obtained, by using the covariant nonlinear dynamical equations for the gravitational and matter fields on the brane, 
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by Maartens • The gauge-invariant formalism for perturbations in the brane world has been developed in |42l Elij . 
The equations governing the bulk perturbations in the case of a general warped Universe have been computed by 
Langlois |44l l45T |. A gauge invariant formalism for metric perturbations in five-dimensional brane world theories, 
which also applies to models originating from heterotic M-theory has been obtained in ^]. Koyama and Soda 01 
obtained a formalism for solving the coupled dynamics of the cosmological perturbations in the brane world and 
of the gravitational waves in the AdS bulk. A closed system of the perturbation equations on the brane, which is 
valid on a large scale and may be solved without solving for the bulk perturbations has been proposed in |4Sj and 
|49| . The perturbations of the brane worlds in conformally Minkowskian coordinates, which enable to disentangle the 
contributions of the bulk gravitons and of the motion of the brane, have been considered in [Hcf. Dorca and van de 
Bruck |5l| proposed a new gauge, in which the full five-dimensional problem of the perturbations is solvable. The 
second order perturbations of the gravitational field induced on the 3- brane have been analyzed by Kudoh and Tanaka 
|52l l53| . The equations of motion for met ric p erturbations in the bulk and matter perturbations on the brane have 
been presented, in an arbitrary gauge, in |54j . The evolution of density perturbations in brane world cosmological 
models with a bulk scalar field has been considered by Brax, van de Bruck and Davis pjfjj . The 1 + 3-covariant 
approach to cosmological perturbations in the brane-world models, and its application to CMB anisotropics, have 
been reviewed recently by Maartens [56j. 

In a previous paper |57| , we studied first order rotational perturbations of homogeneous and isotropic FRW brane 
world cosmological models. Assuming that the rotation is slow, and by keeping only the first order rotational terms in 
the field equations, a rotation equation describing the space dependence and time evolution of the metric perturbations 
is obtained. However, in our previous consideration in the bulk effects related to the tensorial part of the dark 
energy (coming from the 5-dimensional Weyl tensor) have been neglected, and only the role played by the scalar part 
(via the unperturbed field equations) has been considered. The conservation of the angular momentum, following 
from the general energy-momentum conservation equation on the brane, has also not been included in the formalism. 
Moreover, the angular dependence of the angular velocity of the matter on the brane has also been overlooked. 

It is the purpose of the present paper to consider the rotational perturbations of slowly rotating brane worlds, 
with azimuthal symmetry, by taking into account both the tensorial and scalar perturbations in the bulk, produced 
due to the matter perturbations on the brane. We assume that the background geometry of the unperturbed system 
is homogeneous and isotropic, of FRW type. Since the quadratic corrections terms are related, via a consistency 
condition, to the dark energy term, the contributions of the tensorial and scalar parts of the perturbation can be 
consistently included in the formalism. The conservation of the angular momentum of the matter on the brane gives 
a basic relation between the angular velocity uj of the matter and the metric rotation function SI, in terms of the 
equation of state of the matter on the brane and the scale factor of the expanding Universe. In the case of the perfect 
dragging, when lu = f2, the decay of the rotational perturbations is independent of the equation of state of the matter 
and in the large time limit rotation always vanishes. However, the behavior of perturbations is very different for the 
general case u> ^ f2. In such case the analysis of the long time behavior of the angular velocity for several proposed 
physical models of the high-density cosmological fluid, corresponding to different equations of state of the matter 
(Zeldovich causal model or Chaplygin gas) shows that in these cases the rotational perturbations increase in time. 
However, for the Hagedorn equation of state of dense matter, the late time rotational perturbations tend, in the large 
time limit, to zero. For other, more conventional equations of state, like the radiation fluid, pressureless dust, quark 
matter obeying the bag model equation of state and the Boltzmann gas, the rotational perturbations are wiped out 
by the expansion of the Universe. It worth to note that the results for the behavior of perturbations are also specific 
to standard general relativity. 

The present paper is organized as follows. The field equations for a slowly rotating brane-world are written down, 
and the basic rotation equation is derived in Section II. The case of the perfect dragging is discussed in Section 
III. In Section IV the behavior of rotational perturbations is analyzed for several relevant equations of state of the 
cosmological fluid. We conclude our results in Section V. 

II. SLOWLY ROTATING BRANE UNIVERSES 

In the 5D space-time the brane-world is located at Y(X I ) = 0, where X 1 , I — 0,1,2,3,4, are 5-dimensional 
coordinates. The effective action in five dimensions is [30| 

S = J d 5 X^g~ 5 (J^R 5 - A 5 ) + jf ^ d 4 xV=g~ (j^K* - X + L mattor ) , (1) 

with fc| = SttG^ the 5-dimensional gravitational coupling constant and where a; M , fi — 0,1,2,3, are the induced 4- 
dimensional brane world coordinates. R% is the 5D intrinsic curvature in the bulk and is the extrinsic curvature 
on either side of the brane. 
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On the 5-dimensional space-time (the bulk), with the negative vacuum energy A5 source of the gravitational field 
the Einstein field equations are given by 



Gjj = fc§T„, Tjj = -A 59IJ + S(Y) [-X 9lJ + T}J Uer ] , 



(2) 



In this space-time a brane is a fixed point of the Z2 symmetry. In the following capital Latin indices run in the range 
0, 4 while Greek indices take the values 0, 3. 

Assuming a metric of the form ds 2 = (ri/rij + gi^dx 1 dx , with njdx 1 = d\ the unit normal to the \ — constant 
hypersurfaces and gu the induced metric on \ = constant hypersurfaces, the effective four-dimensional gravitational 
equations on the brane take the form [28| 29] : 



G 



where 



1 1 



24 



v (3T at3 T aP - T 2 



(3) 



(4) 



and A = fc 2 (A 5 + fcf A 2 /6)/2, k\ = fefA/6 and Ejj = CiAjBn A n B . Ciajb is the 5-dimensional Weyl tensor in the bulk 
and A is the vacuum energy on the brane. T^ v is the matter energy-momentum tensor on the brane and T = ^ is 
the trace of the energy-momentum tensor. 

For any matter fields (scalar field, perfect fluids, kinetic gases, dissipative fluids etc.) the general form of the brane 
energy-momentum tensor can be covariantly given as 



(p + p)u^u v + ph^ u + 7r^„ + 2q (p u u ) . 



(5) 



The decomposition is irreducible for any chosen 4- velocity u^. Here p and p are the energy density and isotropic 
pressure, tt^ u are the anisotropic stresses on the brane, induced for example by the dissipative properties of the fluid 
(shear viscosity) and g M is the heat flux. h^ v — g^ v + u^u v projects orthogonal to vf". The heat flux obeys the 
condition = g< M >, while the anisotropic stress obeys 7r M „ = 7r <M „>, where angular brackets denote the projected, 
symmetric and trace-free part: 



V< 



w. 



(6) 



The symmetry properties of imply that generally we can irreducibly decompose it with respect to a chosen 
4- velocity field as 



k A 

p - - 5 

^ ~ 1.4 



-UiiUpUv + g^v) + V^v + 2Q {tl u y) 



(7) 



where U is a scalar, Qu a spatial vector and V^v a spatial, symmetric and trace-free tensor. For a FRW model 



Qu = and Vuv = 



Hence the only non-zero contribution from the 5-dimensional Weyl tensor from the bulk 



is given by the scalar term U. The "dark energy" term, E^„, is a pure bulk effect, therefore we cannot determine 
its expression without solving the complete system of field equations in 5 dimensions. However, its expression in the 
bulk is constrained by the relation 



(8) 



The Einstein equation in the bulk also implies the conservation of the energy momentum tensor of the matter on the 
brane, 



T^-v | x =o— 0. 

The rotationally perturbed metric can be expressed in terms of the usual coordinates in the form 

dr 2 o 1 ,„o 



ds 2 - = -dr 



! {t) 



1 — kr 2 



sin 



- 2Cl(t, r, 9)a 2 (t)r 2 sin 2 9 dtdtp, 



(9) 



(10) 



where f2(t, r, 9) is the metric rotation function. Although £1 plays a role in the "dragging" of local inertial frames, it is 
not the angular velocity of these frames, except for the special case when it coincides with the angular velocity of the 
matter fields, k — 1 corresponds to closed Universes, with < r < 1. k = — 1 corresponds to open Universes, while 
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the case k = describes a flat geometry, where the range of r is < r < oo. In all models, the time- like variable t 
ranges from to oo. 

For the matter energy-momentum tensor on the brane we restrict our analysis to the case of the perfect fluid 
energy-momentum tensor, 

= (p + p)u^u v + pgV". (11) 

The components of the four- velocity vector are u M = (1,0, 0, cj) and uj(t,r,6) — dip/dt is the angular velocity of 
the matter distribution. Consequently, for the rotating brane the energy-momentum tensor has a supplementary 
component 

T 03 = {[Q(t, r, 9) - u(t, r, 9)] p - w(t, r, 9)p} r 2 a 2 (t) sin 2 9. (12) 

As a first step to consider rotational perturbations on the brane, we should find a consistent assumption for E^ u 
"induced" by the small rotation. The reason that the small rotation changes the expression of E^ v can be understood 
as follows. In the brane world scenario an Zi symmetry is imposed. Once we introduce a small rotation effect on the 
brane, the bulk geometry should have a "self-tuning" process in order to preserve the Z2 symmetry. Therefore, some 
nontrivia! components of E^ v should be turned on to the order of perturbation. Hence, the first constraint we need 
to check for a consistent expression of E^ is the equation (JHJ). 

For the perfect fluid source, two nontrivial equations of the energy-momentum conservation ©, corresponding to 
the t and <fi components, are 

p + 3(p + p)- = 0, (13) 
a 

fit [{p + p)(fl-u;)a 5 ] = 0. (14) 
Imposing the above equations, the right hand side of © has only one nonzero component as 

D^ = ~d t [pip + p^n-^a 5 ]^ 5 . (15) 

Hence, the constraint JHJ shows that the E^ v should satisfy the following equations 

D li E'* = D„Ei»-=D ll EP a =0, = -^-d t [p(p + p)(n~Lu)a 5 ]a- 5 . (16) 

The last equation simply just shows what we have already mentioned, that is the rotational perturbation does induce 
an effect on the bulk for consistency. Otherwise, the only possible perturbation is the so-called perfect dragging case, 
with Q = lu. 

A possible solution of the constraint (|16fl is E^ v = —(k^/kf)'P^, which, for the perfect fluid energy-momentum 
tensor gives only one non- vanishing component 

k 4 

pt<t> — -p4>t —p^ 0(e), DpEr* = --|a t (Pa 5 )a" 5 . (17) 

Hence the tensor component V is of order of perturbation such that the quadratic terms, by itself or together with u 
or il, can be neglected. Under this assumption, the solution for the equation (jlfijl can be easily obtained 

v = } fp(p+ P m-u J ), as) 

which, as expected, is of the first order of perturbation. Straightforwardly the first order of i?-term is 

hi 

E t4> = E# = -fa 2 r 2 sin 2 9 p{p + p){Sl - to), (19) 

which is the consistently "induced" effect by the rotational perturbation in the brane world model. 

Generally, the above result can be extended to the case including the scalar part of the ^-correction, namely the U 
term contribution. For the case, the first order covariant derivatives are 

DuJ5"* = -H ( d t U + , D^ ir = D,E^ = 0, 



D„E^ = -%d t 



&4 



-W(0-w)o 5 +Pa E 



a -5 + OD M -E^*. (20) 
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Then the constraint (|16|) gives 



U = U§a 



V = 



4 ki , 



(n-w), 



(21) 



Bit, 



K 5 



E b 



Et. 



E, 



4>t 



2 2*2 

a r sm i 



igwn + ^+p)(n- w ) 



(22) 



We assume that rotation is sufficiently slow, so that deviations from spherical symmetry can be neglected. Then to 
first order in f2 the gravitational and field equations © can be decomposed, by neglecting the V^v contribution on 
the background geometry, into the following components 



a 2 3k 



2a a 2 



A - kip - ^p(p + 2p) - 



1 

--U 

3 fc^ 



(23) 
(24) 



3 dc»f7 a 2 o 

a <9r dtdr 



3 a¥ + me ~ °' 



(25) 



2 ^ - f 4 -5*, 



on 



1 <9 2 3cot(9<9ft ~ /2a a 2 
2a 2 fi 



- 2 96» 2 



<96» 



2a 2 fiA + 2A:4a 2 [(fi - uj)p - cop] + yo f{f, > Q, - 2pptu - 2p 2 u>) - 2a 2 
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p(p + p)(n - u) 



(26) 



Equations (|25|l follows from the (r<f>) and (00) components of the field equations and (|26[) is from the (i0)-component. 
Moreover, the last equation can be simplified by applying l|24|) to be 



H , 2 .s 2 (] /4 c7 \ an i a 2 r> 3cot0dft ol2 2 , wo . n 

)^+ (r- 5 ^) — + + — — ^-2fc 4 2 a 2 (p+p)(fi-^=0. 



dr r 2 d6 2 



36 



(27) 



In this equation the correction from the bulk effect exactly cancelled. Hence the form of this equation is the same as 
the version from the standard general relativity. 

The field equations H23|) - (|27|l must be solved together with the conservation equations H13|) - (|14fl . If t ne equation of 
state of the cosmological matter, p = pip), is known, then Eq. I|13f) gives the time evolution of the energy density, 
p = p(a) . Once p is known, the general solution of the background (unperturbed) gravitational field equations on the 
brane can be obtained in the form 



t-tn 



A 2 k 2 
—a , 
3 3 



, , 2 "•" 2/ \ 2 i k§ Uq 

p{a)a + —p (a)a ' 



Z-4 

36' 



3kj a 2 



-1/2 



da, 



where to is a constant of integration. 
The solution of equations (|25H is 

n(t,r,6) = A(r,0)a~ 3 {t), 

where A(r, 9) is an arbitrary integration function. Then the perturbed equations we need to solve are 



(1 - kr 2 ) 



d 2 A 
dr 2 



)kr 



OA 1 d 2 A 3cot6dA 



Or 



de 2 



m 



2k\a 2 {p + p){A- ua 3 ) = 0, 



and 



d t [(p+p)(A-Lua 3 )a 2 ] =0. 



The general solution of Eq. I|31|) is 



A - ioa 3 



F(r,9) 
(p + p)a 21 



(28) 



(29) 



(30) 



(31) 



(32) 
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where F(r, 9) is an arbitrary integration function. With the use of Eq. I|29fl we obtain the following general relation 
between the metric perturbation function f2 and the angular velocity of the matter on the brane: 

W (t, r, 9) = 0(t, r, 9) ^—^ = - j-^. (33) 

Finally, with the use of Eq. (|32|l and of the unperturbed field equations, Eq. I|30[l reduces to 

a - + (i - 5*,) 2^2 + 12!4M + ^M^l _ M|FM) . o. (34) 



III. THE CASE OF THE PERFECT DRAGGING 



As a first case for the study of the rotational perturbations in the brane world we consider the case of the perfect 
dragging, corresponding to the choice F(r, 9) — of the arbitrary integration function. Therefore the angular velocity 
of the matter equals the metric perturbation function and is given by 

w{t,r,9) = A{r,9)a~ 3 (t). (35) 

For an expanding Universe in the large time limit the rotational perturbations always decay, u> tending to zero for 
t — > oo. 

In order to hnd the radial and angular dependence of the angular velocity, we consider two possible functional forms 
of the function A(r, 9). As a first case we take A(r, 9) = B(r)C(9). Then Eq. I|34|) can be separated into the following 
two independent equations 

rPtf dC 

— + 3 C otO- + nC = 0, (36) 

(1 _ kr 2 )r 2 ^- +(t- 5kr) r 2 ^- - nB = 0, (37) 
dr 1 \r J dr 

where n is a separation constant. By introducing a new variable x = cos 9, Eq. Il.''>(:ill becomes 

rp-C rlC 

(l_ x 2 ) l^_ 4a .^l +nC , = (38) 
dx dx 

Therefore the general solution of Eq. I|36[) is given by 

C{9) = d 2 Fi(oi, h; ci; cos 2 9) + C 2 cos 9 2 F 1 (a 2 , b 2 -c 2l cos 2 9), (39) 



where o x = (3 - V9 + 4n)/4, b x = (3 + V9 + 4n)/4, ci = 1/2, a 2 = (5 - y/9 + An)/ A, b 2 = (5 + V9 + 4n)/4, c 2 = 3/2 
and 2^i(a, 6; c; x) = X^fcLo ^°(c)^ fc TT * s tne hypergeometric function [^j. Hereafter in this section the parameters C\ 
and C2 are used to label the constants of integration in the expressions of exact solution. 

In order to solve Eq. I|37[) we introduce first a new variable r\ = r 2 . Hence the equation is transformed into 



,d 2 B „ , „, N dB 

Arj"{l — kn 

with the general solution given by 



4?7 2 (1- krj)^r + 2n(5 ~ 6krj)^f- ~ nB = 0, (40) 
dn z dr\ 



B(r) = C 1 r a 2 F 1 (ki,h;m 1 ;kr 2 ) + C 2 r 2 Fi{k 2 ,l 2 ;m 2 ;kr 2 ), for k = ±1, (41) 



w here a = (-3 - V9 + 4n )/2, /? = (-3 + y/9 + 4n)/2, fc i = -3/4 - V 9 + 4n/ 4, h = 5/4 - V9 + 4n/4, mi = 1 - 
V9 + 4n/2, fc 2 = -3/4 + V9 + in/4, Z 2 = 5/4 + V9 + 4n/4, m 2 = 1 + v / 9T4n/2. For fc = 0, the general solution of 
Eq. (|3"TI) is given by 

B(r) = dr" + C 2 r , for k = 0. (42) 

A second class of solutions can be obtained by assuming that the function A(r, 9) can be represented as A(r, 9) = 
M(r) + N(9). Then Eq. (|34|l yields again two independent equations, 

d 2 N „ dN 

^ + 3cot0- = -m, (43) 

(l-fcrV^+f--5fcrV 2 ^ = m, (44) 
or 2 V r I dr 
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where m is a separation constant. By introducing a new variables N' = dN/d9, Eq. i|43|) is transformed into 

dN' 



dS 



3 cot ON' 



-m. 



Taking integration, Eq. I|45[) gives 



»r/ dN 

N= w =m 



cos 9 1 O 

— — - g cof 5 



sin 



Cr 



sin 



After further integration, we obtain the general solution of Eq. (|43|) in the form 



N{9) = — ( In sin ( 



1 



1 — cos 9 2 cos t 



1 + cos ( 



sin 2 9 



C 2 



The general solution of Eq. I|44|> , corresponding to the three different background geometries is given by 
M(r) = 



l 

1 

6r 3 



mr + \/T + r 2 (2r 2 — l)(msinh 1 r + Ci) — 2r 3 (mlnr + C2)] , 
mr + Vl - r 2 (2r 2 + l)(r7isin _1 r + C\) - 2r 3 (ralnr + C 2 )] , 



d + ^ + f lnr, 



for k = — 1, 
for fc = 1, 
for = 0. 



(45) 



(46) 



(47) 



(48) 



Hence the general solution of the gravitational field equations for a slowly rotating brane with perfect dragging can 
be obtained in a closed form. 



IV. GENERAL TIME DEPENDENCE OF THE ANGULAR VELOCITY FOR u Q 

Generally, the angular velocity of the matter on the slowly rotating brane consists of two time-dependent terms, 
ui(t,r,6) — uji(t,r 7 9) + ui2(t,r,6). The first term, ui(t,r, 6) — A(r,9)a~ 3 (t) always tends to zero for an expanding 
Universe, with an increasing scale factor, so that in the large time limit u>i — > 0. However, the time evolution of the 
second term, u>2(t, r, 9) = —F(r, 9)/(p + p)a 5 essentially depends on the equation of state and the dynamic behavior 
of the background cosmological fluid and geometry. In the following we shall investigate the time dependence of this 
term for a number types of equation of state which could be relevant for the description of the ultra-high density 
matter of which the Universe consisted in its very early stages. 



A. Zeldovich Stiff Fluid 



One of the most common equations of state, which have extensively been used to study the properties of the early 
Universe is the linear barotropic equation of state, with p = (7 — l)p, with 7 = constant G [1, 2]. For this equation of 
state the conservation equation l|13|l can be immediately integrated to give p = poa~ 3j , po — constant > 0. Therefore 
for the second term 0J2(t, r, 9) of the angular velocity we obtain the following time dependence: 

U2 ( t ,r,9) = -^±a 3 '>- 5 . (49) 

IPO 

In order to have decaying rotational perturbations the condition 7 < 5/3 must be satisfied. It is satisfied for the 
case of the pressureless dust, with 7 = 1, and also for the radiation fluid, having 7 = 4/3. But this condition in not 
satisfied for a very important case of the so-called causal limit of the linear barotropic equation of state, corresponding 
to 7 — 2, or the Zeldovich stiff fluid equation of state p = p. For the choice 7 = 2 we obtain u>2(t, r, 9) — —F(r, 9)a/2po, 
showing that an initial rotational perturbation in the ultra-high density cosmological fluid do not decay and the angular 
velocity of the matter is linearly increasing with the scale factor of the expanding Universe. 

The Zeldovich equation of state, valid for densities significantly higher than nuclear densities, p > 10p nuc , with 
Pnuc = 10 14 g/cm 3 can be obtained by constructing a relativistic Lagrangian that allows bare nucleons to interact 
attractively via scalar meson exchange and repulsively via the exchange of a more massive vector meson |60| . In the 
non-relativistic limit both the quantum and classical theories yield Yukawa-type potentials. At the highest densities 
the vector meson exchange dominates and by using a mean field approximation one can show that in the extreme limit 
of infinite densities the pressure tends to the energy density, p — > p. In this limit the sound speed c s = ^Jdp/dp — > 1, 
and hence this equation of state satisfies the causality condition, with the speed of sound less than the speed of light 
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The Zeldovich equation of state can also describe a non-interacting scalar field, with zero potential, with p^ = 
P<6 = 2 /2- Hence the rotational perturbations in a potential-free held do not decay. For a self-interacting field with 
potential U(<j>), the energy density and pressure of the field are p^ — <\> 2 j2 + U (4>) and p^ = (j) 2 /2 — U(<fr), respectively. 
The conservation equation (|13[) becomes in this case 

4> + 3<P- + U' (<P) = 0, (50) 
a 

and has the first integral 

• K- f a 3 U'{cj))dt 
with K a constant of integration. Therefore we obtain for 102 



(51) 



UJ 2 = -F(r, 0) 5 ■ (52) 

K 2 -2Kj a 3 U'{4>)dt + [J a 3 U'{(t>)dt] 

Hence if the potential U(<p) is so that the quantity a{K 2 — 2K J a 3 U' (4>)dt+ [J a 3 U' (fydt] 2 } -1 is a decreasing function 
of time, the initial rotational perturbation is rapidly decaying. 

B. Hagedorn Fluid 

An alternative approach to the equation of state at ultra- high densities is based on the assumption that a whole host 
of baryonic resonant states arise at high density |60| . In the Hagedorn model the baryon resonance mass spectrum is 
given by N(m)dm ~ m a exp(m/m )dm, where N(m)dm is the number of resonances between mass m and m + dm. 
The existing data on baryon resonances show that mo = 160 MeV and —7/2 < a < —5/2 For asymptotically 
large densities the particle number n = f^" N(m)dm ~ mo^exp(p„/mo), where p n is the chemical potential of the 
nuclear matter. The density of the matter is p ~ np n ~ mo/x° +1 exp(yU ra /mo). The pressure can be obtained from 
p = n 2 d(p/n)/dn and is given by the Hagedorn equation of state, 

P=^, (53) 

Po 

where p = 2.5 x 10 12 g/cm 3 

The velocity of sound in this type of matter is c s = [ln(p/po)] _1 ^ 2 [l — (l n (p//°o)) -1 ] ^ |60llolj- For the Hagedorn 
equation of state the speed of sound has the property c s — > for p/po — * 00, in striking contrast with the mean field 
theory approach in which c s — ► 1 . The gravitational collapse of a high-density null charged matter fluid, satisfying the 
Hagedorn equation of state in the framework of the Vaidya geometry was considered in [62]]. A collapsing Hagedorn 
fluid could end either as a black hole, or as a naked singularity. The collapse of Hagedorn fluid to a naked singularity 
is also a possible source of gamma-ray bursts |62| . 

In order to find the time behavior of the rotational perturbations in a Hagedorn fluid, with the equation of state 
given by l|53|l . we have to obtain first the scale factor dependence of the density. For the Hagedorn equation of state 
the conservation equation of the matter on the brane l|13|) takes the form 



p + 3p 



I + (ln-^- 

Pa 



-1 



which can be integrated to give 



Mln£ + 1 



C_ 

n 



= 0, (54) 



(55) 



with C > a constant of integration. For an expanding Universe with an increasing scale factor the energy density 
of the Hagedorn cosmological fluid is decreasing in time. The time variation of the second component of the angular 
velocity of the matter can be written as a function of the time dependent only density as 



2/3 

' ln^ 

"^ e)= - ^ir P , a s ° /3 ' (56) 



Fir. 9) (j) 

( In + I 
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The variation of the time-dependent part f(p/po) = {p/po) 2 ^ 3 m (p//°o)/ (l n (p//?o) + l) 8 ^ 3 
Fig. 1. 



of u>2 is represented in 




2000 4000 6000 8000 10000 
P/Po 

FIG. 1: Variation, as a function of p/po, of the time dependent part f(p/po) of ui2- 

As one can see from the figure, once the Universe is expanding and its density decreases, the angular velocity of 
the matter is also decreasing, and in the large time limit it tends to zero. Therefore the cosmological behavior of 
the rotational perturbations of the Hagedorn fluid is very different to that of the causal stiff (Zeldovich) cosmological 
matter, with the matter angular velocity increasing due to the cosmological expansion. 



C. Chaplygin Gas 



A form of matter which has recently been invoked to account for the recent supernovae evaluations of the cosmo- 
logical constant is the so-called Chaplygin gas |H3 ■ The Chaplygin gas is a form of what is called "k-essence" , and 
satisfies an equation of state of the form p — — 1/p 64] . This equation of state comes from the Born-lnfeld type 
Lagrangian, J (1 — ^/l — Tf- V d^Ad l/ A)d n x, with A the transverse coordinate of domain wall or a (n-l)-brane. The 
energy-momentum tensor satisfies the Hooke law, and, by taking into account a cosmological constant Ao, it follows 
that the pressure and the energy density satisfy the equation of state (p + X )(p — A ) = — 1 [64|]. In the following we 
shall take Aq = 0. Hence the energy conservation equation for a Chaplygin gas on the brane becomes 



P- 



■ s|,-i « 

p) a 



0, 



with the general solution given by 



Po 

«6 



1. 



with po > a constant of integration. 

Therefore the second component of the angular velocity varies as 



w 2 (t, r, 0) 



F(r,0) VaS + p^ 



(57) 



(58) 



(59) 



In the limit of large times a 6 >> po and it follows that w 2 (t, : 



—F(r 7 8)a/ po, that is, similar to the Zeldovich 



fluid, the angular velocity of the expanding Universe is an increasing function of time. 



D. Classical Boltzmann Gas 



Another important example of a rotationally perturbed cosmological fluid we consider next is the classical Boltz- 
mann gas filled Universe. If the cosmological fluid is a collision-dominated classical gas in equilibrium, then the 
thermodynamic parameters of the gas are given by |flfij , [6fiT | 



p = nmj3 



nm = A K 2 {(3)/P, p = A Q [^K^fJ) + 3/3- 2 K 2 (/3)] , 



(60) 
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where n is the particle number density, j3 = mc 2 /ksT ', with ks the Boltzmann constant, T the temperature and 
A = m A g /2ir 2 h? , with g the spin weight of the fluid particles. K n are the modified Bessel functions of the second 
kind. 

With the use of Eqs. (|6U[) the energy conservation equation of the cosmological fluid can be written as 



/3-g(f3)-=0, 



(61) 



where 



6f3[(3Kx(0) + 4K 2 (p)] 



f3 2 K (f3) + SpKiifl) + (12 + f3 2 )K 2 (p) + 3f3K 3 (f3) ' 



In order to obtain Eq. (fr!2"|) we have used the recurrence relation K' n {x) = — [K n -x{x) + K n +i(x)]/2 [&£ 
the scale factor can be expressed as a function of the temperature as 



a(/?)=a exp(/ J) 



(62) 
Therefore 

(63) 



with no > a constant of integration. 

The second component of the angular velocity of the matter on the brane is given by 

(R tt\ F(r,g) ^expf-SjVHl)^) F(r,6) 
where we denoted h{/3) = /3 2 exp [-5 / g _1 (/3) d0\ /\fiKi(J3) + 41f 2 (/?)]. The variation of h(/3) is presented in Fig. 2 



(64) 



4000 
3000- 
^ 2000- 
1000- 

0: 



0.05 0.1 0.15 0.2 0.25 0.3 



P 

FIG. 2: Variation, as a function of /3, of the time dependent part h(/3) of u>2- 
The angular velocity of the Boltzmann gas decreases with increasing j3, so that for small temperatures (large 0) 

LU 2 — * 0. 

E. Quark Matter 

The rotational perturbations always decay in the important case of quark matter, satisfying the bag model equation 
of state p = (p — 4_B)/3 |6jJ, where B — 10 14 g/cm 3 is the bag constant 68]. Quark matter is assumed to have played 
an important role in the early evolution of the Universe, and a first order phase transition from quark to hadronic 
matter took place when the temperature was around lOOMeF. It is also possible that inhomogeneities in the baryon 
number density were produced during this transition, and perhaps they persisted even to the time of nucleosynthesis, 
which could alter the abundance of light elements 67] . For the bag model equation of state of the quark matter the 
density-scale factor dependence is given by 



B 



(65) 
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and consequently u>2 becomes 



„,(«,,.,„ —J^iffli. (66, 



Therefore in the large time limit, for a — > oo, the rotational perturbations in the cosmological quark matter tend to 
zero. 

In the general case, due to the dependence of the rotation equation on the arbitrary function F(r, 9), the solution of 
the rotation equation cannot be obtained. The particular choice F(r,9) = A(r,0) leads again to a radial and angular 
distribution of the matter which is very similar to the perfect dragging case. Moreover, if A is represented as a product 
of two r and 6 only dependent functions, A(r, 9) = X(r)Y(9), then the angular velocity w is also a separable function 
of all variables. But if A is the sum of two independent functions, A(r,9) = X(r) + Y(9), the angular velocity is a 
non-separable function. In these cases the functions X(r) and Y{9) can be obtained in terms of the hypergeometric 
function. 



V. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have considered the evolution of the small rotational perturbations in the brane world 
cosmological model. The evolution of the angular velocity of the matter is related to the metric perturbation function 
via an equation derived from the energy-momentum conservation on the brane. This equation is in fact the same in 
both standard general relativity and brane world cosmology. Once the evolution of the background non-perturbed 
geometry is known, the time dynamics of the metric rotation function f2(f , r, 9) and of the angular velocity w(t, r, 9) 
of the matter on the brane is uniquely determined by the unperturbed field equations and, from a physical point 
of view, by the equation of state of the ultra-high density cosmological fluid. The time behaviors of ui(t,r,9) and 
fi(f , r, 9) are independent on their spatial distribution. For to = Q, in the large time limit the angular velocity of the 
matter tends to zero, showing that the initial small rotational perturbations of the brane world are smoothed out, 
due to the expansion of the Universe. This result is also independent of the geometry (flat, open or close) of the 
space-time. However, despite the fact that the presence of the dark matter term U = Uo/a 4 is not essential for the 
rapid decay of the rotational perturbations of the matter on the brane, the effects induced by the five-dimensional 
bulk still contribute to the decay (or amplification) of the first order rotational perturbations via their effect on the 
unperturbed background geometry. 

A very different situation occurs for u> ^ fi. In this case the dependence of the angular velocity on the equation 
of state leads to two distinct classes of behaviors. For some equations of state (radiation fluid or quark matter), the 
rotational perturbations decay and this behavior is consistent with the observational constraint that our present day 
Universe is rotating very slowly. However, some equations of state which try to model the thermodynamical properties 
of the super-dense matter at the very early stages of the evolution of our Universe, like the Zeldovich or Chaplygin 
gas equations of state, do not lead to an observationally consistent description of the early Universe. In particular, 
the causal Zeldovich equation of state with 7 = 2, which has been used, for example, to determine the maximum mass 
of the neutron stars |6(| , does not lead, once applied to a cosmological framework, to a description of the rotational 
perturbations consistent with the observations, that is, it does not satisfy the basic observational requirement of a 
very small rotation in the large time limit. Hence, if still one assumes these equations of state for the early Universe, 
a physical mechanism supplementary to the expansion is also needed to suppress rotation. 

However, if the brane Universe experiences an inflationary period, than any rotational perturbation is reduced 
exponentially to a very low level, the decrease of cosmological vorticity being of the order of 10 -145 01 • Hence, 
any growth of the rotational perturbations during a maximally stiff phase in a post-inflationary period is unlikely to 
increase them back up to observationally significant levels. 

The equation of state of the Hagcdorn fluid is consistent with the requirement of a non-rotating Universe. Therefore 
the study of the perturbations of cosmological models could also provide some insights in the equation of state of the 
ultra-high density nuclear or sub-nuclear matter. The actual very tight limit on the rotation of the Universe imposes 
very strong constraints on the initial equation of state of the cosmological fluid. 

The spatial and angular distribution of the angular velocity depend on two arbitrary integration constants C,*, i = 
1,2. In principle the values of these constants could be determined by fixing the value of A(r, 9) at the center, and 
from the actual value of the angular velocity of the Universe. 
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